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Abstract 

An inequality of the eigenvalues of the reduced density matrix p 2 at finite 
temperature in the Hubbard model is obtained by means of the Bogolyubov 
inequality. The quasi-average of S + in a simple symmetry-breaking pertur- 
bation of the Hamiltonian for a bipartite lattice is shown to be zero. 
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^ Hubbard models and its extensions such as the SO (5) model ]]]] are currently believed to be 
able to account for the high-T c superconductivity of cuprates. The possibility of existence at 
finite temperature of s-wave pairing, generalized //-pairing and d x 2_ y 2 in one and two dimen- 
sional Hubbard model and Anderson model seems to have been be ruled out [^J (of course, 
the conclusions are only valid with respect to the special kind of symmetry-breaking per- 
turbations, there may exist some other kinds of perturbations for which the quasi-average 
of the order parametres do not suffer from those limitations, and thus are possibly non- 
vanishing). On the other hand, there do exist some states in some Hubbard models exhibit 
*On leave of absence from the Physics Department, Shanghai University, 201800, Shanghai, China 



superconductivity M. Therefore, the relevence to superconductivity of Hubbard models 
still calls for further investigations. A convenient concept is off-diagonal long-range order 
(ODLRO) discussed in detail in0 which can be studied by //-pairing 0and its variants. 
The existence of ODLRO can be recognized if the largest eigenvalue of the two-body re- 
duced density matrix p 2 is of the order of the total number of the fermionsf^]. A novel 
property of ^-pairing is that it constutites a pseudo-spin algebra. Lieb used this algebra and 
partial particle-hole transform and exactly proved for the first time that ferro-magnetism 
may exist in itinerant electrons 0. Base on Lieb's results, it is shown that for some specific 
lattice structure and electron-filling, the ground-state of negative-f/ Hubbard may support 
ODLRO |Q. In this letter, we first prove an inequality of the eigenvalues of p2 by means 
of Bogolyubov inequality. Then we show by means of fluctuation- dissipation theorem that 
the grand canonical quasi-average of S + vanishes for both positive and negative-f/ Hubbard 
models. 

The grand-canonical Hamiltonian for the Hubbard model is 

K = E E Uj<Lc j(T + n^nn - fiN (1) 

(ij) ° i 

where c\ a and q ct are the creation and annihlation operators of the electrons with spin 
a =|, I at site i. The hopping matrix {tij} are required to be real and symmetric. The 
number operators are n icr = c\ a Ci a , while the U denotes the on-site Coulomb interaction. It 
is further assumed that the the lattice A is bipartite in the sense that it can be devided into 
sublattices A and B, i.e. A = A U B, such that = o whenever {ij} G A or {ij} G B. As 
usual, the spin S and pseudo-spin S are defined as follows 

S + = E 4 T <*> S ~ = E 4|Q T , S z = \ E(4 T Q T - c\ iCii ) = \{N,- N t ) (2) 
ieA ieA ieA 

S + = E e< S)^ c ih S~ = E e (O c !| c Ir> S* = - E(l - - n n) (3) 

ieA ieA ieA 

where e(i) = lwhen ieA and —1 when i G B .Both the spin and the pseudo-spin operators 
constitute SU(2) algebra and they commute with each other, i.e. [S, S] = , so they form 
an SU(2)® SU(2) algebra. It is not difficult to show that [H, S 2 ] = [if, S] = [H, S M ] = 0. 
Therefore, Hubbard model enjoyes SU(2) g) U(l) ® U(l) symmetry. Yang and Zhang[^] 
showed that ODLRO exists whenever the expectation value of S 2 — Si is of order A 2 , where 
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N\ is the number of the sites of the lattice considered. 

The well-known Bogolyubov inequality U says that for any two operators C and Q, the 
following inequality always holds 



<CC +cc> - (3\<[Q(t),[QKt),K]]> 



2 

(4) 
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where < O > denote the statistical average of the mechanical quantity O. This inequality 
has been utilized in[[|. We employ it here to prove a relation the eigenvalues of the reduced 
density p2 should satisfy Let us recognize some fundamental properties of the eigenvalues 
before presenting the relation. Firstly, it is not difficult to show that the onset of ODLRO 
can be guaranteed if the largest eigenvalue of either < e(j)r)je(i)r)i > or < r/jr/j >, where the 
T) pairing is defined as rji = CjfCj|, is of the order of N\. In fact, we have the following lemma. 
Lemma 1 The matrices < e(J)rjje(i)r)i > and < ry^rji >. have the same eigenvalues. 
Proof Suppose J2i < VjVi > u i = Xuj, then we have 

< e (j'HW(0 > = e(j) J2 < VjVi >Ui = Xe(j)uj (5) 

i i 

and this proves the lemma. Q.E.D. 

Secondly, the eigenvalues are non-negative. 
Lemma 2 The matrix < e(j)r]je(i)r]j > is positive semi- definite. 
Proof For an arbitrary vector m, define O = Ya e iVi u ii then 

J2 u* < ■n]e{j)e{i)r H > u % =< ] >> (6) 

Therefore, all the eigenvalues are non-negative. Q.E.D. 

Beside the above two lemmas, we will use the next lemma in the following. 
Lemma 3 Let A be an N x N matrix whose entries satisfy a mn = a m - n . Then all the 
eigenvalues of A are given by X q = J2m=i ^n=i a mne tq(ja ~ n \ where q = j^k with k = 

o,i,...,jv-i. 

This lemma has been used in and |TD[ and a proof is presented in [HJ . Now we are ready 



to present our first result 

Theorem 1 The eigenvalues of< rjjTji > can be denoted by wave vectors k as Ak, k belong to 
the first Brillouin zone. If the system is translationally invariant and the lattice is bipartite, 
the following inequality holds 



where x = 1 — p e . 

Proof For our purpose, we choose 

Q = S~ (8) 

and 

Ck = £ c^ ( (j)r h = X e^6(i)c jTCji (9) 
Making use of the relation [r)},r)j] = Sij(n^ — 1), we have 

[Q,cy = £e^K-i) (io) 

j 

Since 

[S-,K\ = (2iM-U)S~ (11) 

we have 

[Qt,tf] = (tf_ 2/ i)S + (12) 
[Q,[Q\K]]=2(2fi-U)S z (13) 

By direct calculation, we have 

CtC k + C k c£ = 2 £ e*-^-^)c(<)cO>^ - ]T(n, - 1) (14) 

ij i 

Therefore, using Bogolyubov inequality, we have 

2£e lk - (R <- R ^(*) e (j) < ^ > -£ < m - i > 

ij i 

2 | < Ej e <k,R ^(nj - 1) > | 2 _ 2 Y*j e ik - (R ^ R ^ < ^ - 1 >< rij - 1 > 
" P \2(2fj,-U)\\ <S Z >\ ~ P 2\2fi - U\ ■ \<l-n j >\ 
Since we consider the translationally invariant system, < rii >= N e /N\ = p e , accordingly 

where we have used that e lk '( Rl ~ R ^ = N A 5 kt0 - Or 

2 $>*■<*" »•• < 6(j>t c (o^ >> - ^ ( 17 ) 

On the other hand, lemma 3 tells us that the eigenvalues of the N A xN A matrix < e(j)rjje(i)r]i > 



are 



Ak: AT A - 



1 ^ elk .(R^R j)<e(j)r? t e(?k> (lg) 



A 



and this proves the theorem. Q.E.D. 

We make some remarks here. In the case of at half-filling, x — 0. in this case, the 
theorem provides no new information because lemma 2 has told us that all the eigenvalues 
are non- negative. It should be mentioned that Tian showed in|l(J that at half- filling, for 
the negative— U case, the lowest eigenvalue is zero in the ground-state. Our result is in 
accordance with Tian's, though the zero-temperature limit of a grand-canonical average 
does not necessarily coincide with the ground state expectation value in general such as| 
and 0]. But, if the system is at over half-filling, x will be negative and the r.h.s. of (7) is 
positive, in this case, the theorem tells us that there exists lower positive bound to each 
eigenvalue. Unfortunately, the lower bound to the largest eigenvalue the theorem provides 
is only a constant other than a number proportional to N e . 

Our second result is 

Theorem 2 Assuming the lattice is bipartite, for the U(l) symmetry-breaking Hamiltonian 
H = E E tijclcj* + U E rum + AS + + A*S- (19) 

(ij) <r i 

we have for the grand canonical average of S + 

lim lim — < S + >= (20) 

A^O V A ^oo jV A 

Proof The grand canonical Hamiltonian is K = H — fiN. We use the fluctuation-dissipation 
theorem as in [U\. to calculate < S + S~ >. The evolution equation of the double-time 



Green function C S |S + ^> w . isJT2| 



w<S-|S+> w =< [ST, S+] > + < [S-,A^]|S+> W (21) 
It can be calculated directly that 

[S-,K\ = (2fi-U)S- -2AS Z (22) 

Accordingly, we have 

uj < S-\S + > w = -2 < S z > +(2fi - U) < S-\S + > -2A < S Z \S + > w (23) 
Using again the equation of motion, we can obtain 

u < S-\S + > w = - < S + > -A < S + \S + > w +A* < S-\S + > w (24) 



and 

to < S + \S + >^= -(2/i - U) < > w +2A* < > w (25) 

These three equations eq.(23)-eq.(25) are enough to determine the three unknowns <C S~\S + 3> 
S + \S + ^$> u , and <C ^> w . Especially, we have 

« S-\S+ »„= - + — ^— + (26) 

UJ UJ — UJ 2 tO — Co>3 



where 



Cl = [2A(2a* - U) < ~S + > -4|A| 2 < S z >} AW _\^_ U)2 (27) 



c 2 = — -{2A <S + > (lu 2 + 2fM-U)- 2[lu 2 (lu 2 + 2fM-U) + 2|A| 2 ] < S z >} (28) 

UJ 2 {UJ 2 - UJ 3 ) 

c 3 = — — ^ -{2A < 5+ > (cj 3 + 2(jt-U)- 2[u 3 {uj 3 + 2p-U) + 2|A| 2 ] < S z >} (29) 

^3(^3-^2) 



W2i 3 = ± v /(2 /U -f/)2-4|A|2 (30) 
Thus from the fluctuation-dissipation theorem 



< S + S- >= — du — 1 " W+ Z ^ 1 (31) 

2vr 7-oo - 1 v ; 

we know that C\ must vanish, otherwise, there will be < S + S~ >= 00 for finite N\. Thus 

1 ~ 2A* 1 

TT < S+ >= ~^T-r4- <S Z > (32) 
N A 2/i -UN K K J 

Taking the limit in (20), one can immediately arrive at the conclusion Q.E.D. 

According to Bogolyubov |l3| , symmetry-breaking properties should be discussed by 
quasi-averages, i.e. giving some symmetry-breaking perturbation and letting it vanish finally. 
Theorem 2 states that the quasi-average of S + is zero. This does not mean the pairing 
liniA^o limjv A ^oo < QfCjj >= 0. If the two sub-lattices are homogeneous, denote s\, B =< 
CjfC^ > for i G A/B. Theorem 2 implies that pas\ — Pbs% = 0, where pa/b = Na/b/N\. 

We make some conclusional discussions. It can be understood that we can not obtain a 
lower limit of the order of N\ for the largest eigenvalue of the reduce density p 2 by means 
of the Bogolyubov identity (4) because the times of sum J2i m both the numerator and 
the denominator in the right hand side of (4) are equal . Thus the value will not be of 
order N\. As the conclusions in 0] hold only for the specific quasi-averages therein, our 
theorem 2 refers to only a special way of symmetry breaking, or a special quasi-average too. 



It can not surely rule out the possibility for non- vanishing < S + > for other quasi-averages. 
Similar discussions seem to apply to other models such as Kondo lattice model and periodic 
Anderson model. Finally, we would like to stress that both theorems do not rule out the 
possible existence of ODLRO in Hubbard models. 
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